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Abstract 

We calculate the lowest order QCD amplitude, i.e. the quark exchange contribution, 
to the forward production amplitude of a pair of longitudinally polarized p mesons 
in the scattering of two virtual photons r y*(Qi)l*(Q2) —> P°lP°l- We show that the 
scattering amplitude simultaneously factorizes in two quite different ways: the part 
with transverse photons is described by the QCD factorization formula involving 
the generalized distribution amplitude of two final p mesons, whereas the part with 
longitudinally polarized photons takes the QCD factorized form with the — > p\ 
transition distribution amplitude. Perturbative expressions for these, in general, 
non-perturbative functions are obtained in terms of the p— meson distribution am- 
plitude. 
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1 Introduction 



The exclusive reaction 

l*(qih*(q2)^P° L (h)p L (k 2 ). (1) 

is a beautiful theoretical laboratory for the understanding of factorization properties of high 
energy QCD. In previous studies we emphasized the gluon exchange contribution which 
dominates in the very high energy (small-x) regime, both at the Born order [1] and in the 
resummed BFKL approach [2]. At lower scattering energy one can expect a non negligible 
contribution of quark exchange, since in the strong coupling g expansion, quark exchange 
processes appear at lower order, namely at order g 2 , than the gluon exchange ones. The 
study of these quark exchanges processes is the main motivation of the present work. The 
Born order contribution with quark exchanges is described by the same set of diagrams 
which contribute to the scattering of real photons producing pions, e.g. 77 — ► 7r + 7r~, at large 
momentum transfer, studied long ago by Brodsky and Lepage [3] in the framework of the 
factorized form of exclusive processes at fixed angle [4], in which mesons are described by 
their light-cone distribution amplitudes (DAs). This is illustrated in FigJTJ In this sense our 




Figure 1: The amplitude of the process 7*(Qi)7*(Q2) — ► P^i^i) Pl^) in the collinear factorization. 

present study can be seen as a complement of Ref. [3] for the case of the scattering with 
virtual photons, i.e. with both, transverse and longitudinal polarizations, and in the forward 
kinematics (see also Ref. [5]). The virtualities Q 2 = —q 2 , % = 1,2, supply the hard scale to 
the process (JIJ which justifies the use of the QCD collinear factorization methods and the 
description of p mesons by means of their distribution amplitudes. 

We calculate in Sec. El within this scheme the scattering amplitude of (JTJ) at Born level. 
At this stage we do not impose any additional conditions on the magnitudes of photon 
virtualities Qi. Next, we turn to the study of two particular kinematical regions: (1) the 
region where the squared invariant mass of the two rho's W 2 is much smaller than the 
largest photon virtualities, namely Q\ 3> W 2 (or Q\ ^> W 2 ), with Qi and Q 2 being not 
parametrically close, see Sec. 0] for more details, and (2) the region where photon virtualities 
are strongly ordered, that is Q\ » Q\ (or Q\ » Q\). We will show in Sec. 0] that 
in the region (1) the amplitude with transverse photons factorizes in a hard subprocess 
and a Generalized Distribution Amplitude [6,7], up to corrections of order W 2 /Q\ (resp. 
W 2 /Q\). In the region (2) the amplitude with longitudinal photons factorizes in a hard 
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subprocess and a Transition Distribution Amplitude [8], up to corrections of order Q\jQ\ 
(resp. Q\/Q 2 ), as is shown in Sec. O These two domains have a non empty intersection as 




Figure 2: Different kinematical regions. In the domain (1), denoted with vertical lines, the QCD 
factorization with a GDA is justified. This region is the union of two disconnected domains: the 
lower right one is given by the conditions 1 — Q 2 /s > c(l — Q\/s) and Q 2 /W 2 > c while the 
upper left one is given by 1 — Q 2 /s > c(l — Q\j s) and Q\jW 2 > c, c being arbitrary large. In 
the domain (2), denoted with horizontal lines, the QCD factorization with a TDA is justified. It 
corresponds to Q\jQ\ > c or Q 2 /Q 2 < 1/c. In the intersecting domain, both factorizations are 
valid, while in the region without any lines, no factorization neither in terms of GDA nor TDA 
is established. For illustration, we choose c = 3. The grey band represents the domain where the 
Born order amplitudes with transverse and longitudinal photon, calculated directly as in Sec. El have 
comparable magnitudes. In the upper (lower) corner the transverse (resp. longitudinal) polarizations 
give dominant contribution. Below the dashed line is the perturbative Regge domain. 

shown in Figj2j In this intersection, we get two factorisation formulas for both polarizations 
of the photons, 7^7t — > PlPl and 7272 — > PlPl- Fig. [2] illustrates also that the collinear QCD 
factorization with GDA or TDA is not demonstrated by our analysis in the limited region 
where the virtualities Q 2 are parametrically close. In this figure, the different domains where 
the transverse or longitudinal amplitudes dominates are displayed. Finally, Fig. [21 shows the 
perturbative Regge domain, corresponding to the large W 2 limit. 

In this paper we concentrate ourselves, for simplicity, on the case of longitudinally polar- 
ized rho mesons. In this spirit, the longitudinal rho pairs and the pion pairs lead to similar 
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results. Our study is closely related to the perturbative limit of the two-pion light-cone 
distribution amplitude performed in Ref. [5] for 7^7t collisions in the GDA limit. 

Thus, the question we investigate here is how the known QCD factorization formulas 
emerge from a Born order calculation in both the longitudinally polarized and transversally 
polarized photon cases. As a by product, we get a perturbative expression for the non- 
perturbative hadronic objects - the p p GDA and the 7* — > p TDA - in terms of the p— meson 
distribution amplitude. 

2 Kinematics 

Let us fix the kinematics, which is quite simple in the forward regime that we are investigating. 
The inclusion of transverse momentum would not change dramatically our result but would 
necessitate the introduction of more complicated tensorial structures. We use a Sudakov 
decomposition with two lightlike vectors p\ and p 2 with 2p 1 .p 2 = s and write the photon 
momenta as 92 

Ql Q2 /r>\ 

ft = Pi P2 Q2=P2 Pi , (2) 

s s 

and the final meson momenta as 

h = (1 - ^) Pl k 2 = (1 - ^)p 2 . (3) 

The positivity of energy of produced p's requires that s > Q\. The usual invariant W 2 is 
defined as 

W 2 = (ft + q 2 f = (h + k 2 f = 8 (1 - - 5!) , ( 4 ) 



s 



or 



s = \[Q\ + Ql + W 2 + \{Ql Ql -W 2 )] , (5) 

with X(x, y, z) = \J x 2 + y 2 + z 2 — 2xy — 2xz — 2zy , while the minimum squared momentum 
transfer is 

t mm = (q 2 - hf = ( gi + h) 2 = -9M . (6) 

Note that, contrarily to the case studied in Ref. [3,5], t m i n may not be large with respect to 
Aq CD . depending on the respective values of Ql, Q\ and W 2 . 



3 The Born order amplitude 

The scattering amplitude A of the process ((H) can be written in the form 

A = T""eMe v (q 2 ), (7) 

where the tensor T^ v has in the above kinematics a simple decomposition which is consistent 
with Lorentz covariance and electromagnetic gauge invariance 

T» v = \& v (rV^) + (pf + %" 2 W 2 + —Pi) - 2 (T^ P2aPip ) , (8) 
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and where gt^ v — g^ v — (piP 2 + PiP 2 ) I (pi-Pi)- The first term on the rhs of Eq. JHJ) contributes 
in the case of transversely polarized photons, the second one for longitudinally polarized 
virtual photons. The longitudinally polarized p°— meson DA (p(z) is defined by the non-local 
correlator (for simplicity of notation we omit the Wilson line) 

1 

(p° L (k)\q(x)Yq(0)\0) = 4 & J dz*W»H{z) , for q = u, d , (9) 



proportional to the coupling constant f p . 

The Born order contribution to the amplitude (jZj is calculated in a similar way as in the 
classical work of Brodsky-Lepage [3] but in very different kinematics. In our case the virtual- 
ities of photons supply the hard scale, and not the transverse momentum transfer. Moreover, 
we consider the forward kinematics for simplicity. Relying on the collinear approximation, 
the momenta of the quarks and antiquarks which constitute the rho mesons can be written 
as 

£ 1 ~ zik 1: £ 2 ~ z 2 k 2 

£1 ~ z x k u £ 2 ~ z 2 k 2 . (10) 

The number of possible diagrams is 20. They can be organized into two classes. The first 
class corresponds to the diagrams where the two virtual photons couple to two different quark 
lines. In this way, one can build 8 different diagrams, as illustrated in Fig03 The second 




Figure 3: Feynman diagrams contributing to Mh, in which the virtual photons couple to different 
quark lines. 
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Figure 4: Feynman diagrams contributing to Mjj, in which the virtual photons couple to a single 
quark line. 

class of diagrams corresponds to the one where the two virtual photons are coupled to the 
same quark line, resulting into 12 different contributions, as illustrated in FigQ] 

Let us first focus on the case of longitudinally polarized photon. Their polarization read 

e ll = TT^ 1 + P 2 and e \\(lV = TT^ 2 + Pi- ( n ) 

Vi s y^>2 s 

Due to the gauge invariance of the total amplitude, when computing each Feynman diagrams, 
one can make use of the fact that the longitudinal polarization of photon 1 (resp. 2) can be 
effectively considered as proportional to P2 (resp. pi), in accordance with the structure of 
the second term of Eq.(JEJ). One then readily sees that in the kinematics we are investigating, 
the only diagrams which are non zero in the Feynman gauge are the diagrams (16) <8> (2a), 
(3a) ® (46) on the one hand, and (s2), (s2') on the other hand. Thus, only the four diagrams 
illustrated in FigJHl contribute. 

In the case of transversally polarized virtual photons, described by the first term in Eq.(j§|). 
one needs to contract the two polarization indices through the two dimensional identity tensor 
g^ v . The four non vanishing graphs corresponding to the coupling of the photons to different 
quark lines are (la)<8>(26), (16)(g>(2a), (3a)®(46) and (36)®(4a). When considering the graphs 
corresponding to the coupling of the photons to the same quark line, the above mentioned 
contraction kills the graphs (s2), (s2') and (s5), (s5'), that is the graphs where the gluon is 
emitted from the quark connecting the two virtual photons. In this second class of diagrams, 
8 diagrams thus remain. The whole set of twelve diagrams to be computed is shown in FigEl 
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Figure 5: Feynman diagrams contributing to Mh in the case of longitudinally polarized virtual 
photons. 



It should be noted that in the peculiar kinematics which is chosen here, the only diagrams 
which contribute both to the longitudinal and to the transverse virtual photon cases are 
(16) ® (2a), (3a) ® (46), that is the two first diagrams of FigJSland FigG3 

The scalar components of the scattering amplitude (JHJ read : 



ia(3 



9Ta(3 



e 2 (Ql + Q 2 d )9 2 C F p p 



x 



2 1 



Q_ 

s 



4iV cS 

_Ql 

s 



dzi dz 2 <j)(zx) 4>(z 2 ) 



{Z2 + Z2^) 2 (Z 1 +ZS) 2 



+ 



Z 2 Z\ Z\ Z2 



1 



1 



{Z 2 + Z 2 - 
1 



} h 2 (zi + zS) 2 



(12) 



+ 



1 - — \z 2 + z 2 — 

s 2 f 2 p C F e*g\Ql + Ql) 
8N c QjQl 

(1 -¥)(!-¥) 



Z2 + Z 2 ^. 



s 



dz\ dz 2 (j)(zx) 4>(z 2 ) 



(13) 



+ 



1 1 
+ — — + — 



(Zi+Z^iZz+Z^) (z 1 +z 1 ^)(z 2 + z 2 ^) *2*1 z lZ2 

where Q u = 2/3 (Q d = -1/3) denote the charge of the quark u (d), C F = (iV c 2 - l)/(2N c ) 
and N c = 3. 

The above results are obtained for arbitrary values of the photon virtualities Qi. A closer 
look into formulas ([7l lT3|) leads to the conclusion that all integrals over quarks momentum 
fractions Zi are convergent due to the non-zero values of Qi. Further inspection of the 
amplitudes (f71ITTS|l reveals that the transverse photon part has a behaviour like 1/W 2 while 
the longitudinal photon part behaves like l/(QiQ 2 ). This shows that their relative magnitude 
is different in each region (1) and (2) discussed in Sec. 1. 
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(si) (s3) (s4) (s6) 




(si') (*3') (s4') (s6') 



Figure 6: Feynman diagrams contributing to Mh in the case of transversally polarized virtual 
photons. 
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ItIt 



p° L Pi in the generalized Bjorken limit 



Let us first focus on the region where the scattering energy W is small in comparison with 
the highest photon virtuality Qi 



W 2 _s_ 

Ql ~Qi 



Ql 



Q 



l 



Ql 



<c 1 



(14) 



which leads to the kinematical conditions very close to the ones considered in Ref. [6, 7] for 
the description of 7*7 — > nn near threshold. In Ref. [6, 7] it was shown that with initial 
transversally polarized photons, the scattering amplitude factorizes at leading twist as the 
convolution of a perturbatively calculable coefficient function and a generalized distribution 
amplitude (GDA). We will recover a similar type of factorization with a GDA of the expression 
(fT2|) also in the case of our process ((H), as illustrated in Fig0 provided Qi and Q 2 are not 







1 (gn da 



2>h (gr da 
^rv_yp(A; 2 ) 

Figure 7: Factorisation of the amplitude in terms of a GDA. 




parametrically close, i.e. 



1 _ 91 « 1 _ 91 . 

s s 



(15) 



Indeed, in the scaling limit (fT4l) with ((TBI) , the contribution of the four cats-ears diagrams 
(i.e. the first line of FigEJ), which corresponds to the second line of (TP2| . is subdominant and 
the dominant contribution is given by the last term in (fT2|) . i.e. 



e 2 (Ql + Q 2 d )9 2 C F p p 




dzi dz 2 <t*(zi) 4>(z 2 ) 



(16) 



The virtuality Q 2 plays in the expression (fTKl) the role of a regulator of the end-point singu- 
larities (compare with [7]). Before recovering each factor of the factorized equation (fTK|) by 
a direct calculation, let us discuss the physics which is behind. In Fig® the diagrams which 
contribute in the scaling region ([Hj) with (fl~5|) are (s3), (s3') and (s6), (s6'). The sum of 
these four diagrams factorizes into a hard part convoluted with a soft (but still perturbative) 
part. Indeed, the typical virtuality of the hard quark connecting the two virtual photons is 
of the order of Q\ for graphs (s3), (s3') and s for graphs (s6), (s&), while the virtuality of 
the quark connecting the emitted gluon to the virtual photon is of the order of W 2 , which is 
negligible in the scaling region with respect to Q\ and s. 
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Let us now recall the definition of the leading twist GDA for p\ pair. We introduce 
the vector P — k± + k 2 ~ pi, whereas the field coordinates are the ray- vectors along the 
light-cone direction n M = p 2 /(pi-p 2 ). In our kinematics the usual variable ( = (kin)/(Pn) 
characterizing the GDA equals £ ~ 1. Thus we define the GDA of the p° L pair & q (z, £, W 2 ) 
by the formula 



(p° L (h) P° L (k 2 )\q(-a n/2)i 



dy n u A v (y) 



q(a n/2)\0) 



q = u,d. 



(17) 



Now we calculate the GDA Q q (z, £, W 2 ) in the Born order of the perturbation theory. 
First we show that the gluonic Wilson line does not give a contribution in our kinematics. 
For that expand the Wilson line and the S-matrix operator with quark-gluon interaction in 
(|T7|) at the first order in g. We obtain (up to an irrelevant multiplicative factor) 



n,9 2 / d i v{pl{k 1 )pl{k 2 )\T[q{-an/2)^ 



dyn u A v (y) 



q(an/2)q(v)A(v)q{v)}\0) (18) 



After applying the Fierz identity and ordering quark operators into two non-local correlators 
defining DA of p— mesons we obtain 

^(PL(^i)l?(-W2)7 5 gM|0)(p° L (A; 2 )|g(t;)7 a g(W2)|0)Tr[7V7 CT 7 /3 ]n /3 + (h <- k 2 ) , (19) 

where we omitted the gluon propagator coming from the contraction of the two gauge fields. 
The original Wilson line results at this order in the presence of the vector rip in the above 
expression. The only nonvanishing contribution at leading twist takes the form 



(pi(^i)l?(-W2)?ig(«)|0)(pi(A;2)|g(i;)fig(an/2)|0)Tr^i^] + (h <-> k 2 ) 
This is illustrated in FigJHl 



(20) 




+ 




Figure 8: Wilson line contribution to the GDA, expressed as a convolution of a hard part with the 
DAs. 

It equals zero since in our kinematics one of the matrix elements defining DA of p— meson 
vanishes. 
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Figure 9: Non vanishing contributions to the hard part of the GDA at Born order. The quark and 
gluon bold lines correspond to propagators . 



The remaining contributions to the correlator in (fTTj) at order g 2 are illustrated in FigEl 
It leads to the result 

$ pLPL (z,( « 1, W 2 ^ - p 



ft 9 2 C F 



2N r W 2 



dz 2 (f){z) (p(z 2 ) 



1 1 

zz 2 zz 2 



211 



The hard part T H of the amplitude corresponds to the diagrams shown in FigUHl In the case 




Figure 10: Expansion of the hard part Tjj at g 2 order. The bold lines correspond to quark propa- 
gators. 



of a quark of a given flavour it equals 

T H (z) = -4e 2 N c Q 



"' 1 z + z^ z + z^ 

s . 



(22) 



The Eqs. (|2"T1 122J) taken together with the flavour structure of p° permit to write ifTfijl in the 
form 



(Ql + Ql) [ dz 

J \z + z 



?2 _l ~Q~2 



^CwlX) , (23) 



z + z- 



which shows the factorization of T a/3 g Ta/3 into the hard part and the GDA. The Eq. (|23| is 
the limiting case for ( — > 1 of the original equation derived by D. Miiller et al [6]. 



5 ~^ PlPl m ^ ne generalized Bjorken limit 

To analyse the case of longitudinally polarized photon scattering, let us now turn to another 
interesting limiting case where 

Qi » Ql- (24) 
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In this regime, it has been advocated [8] that the amplitude with initial longitudinally po- 
larized photons, should factorize as the convolution of a perturbatively calculable coefficient 
function and a 7 — > p transition distribution amplitude (TDA) defined from the non-local 
quark correlator 



"> Z ^-ixP+2 



2tt 



< p( P 2)\q(-z'/2) 1 + q(z-/2))\ 7 (g 2 ) > , (25) 



which shares many properties (including the QCD evolution equations) with the general- 
ized parton distributions [6, 9] succesfully introduced to describe deeply virtual Compton 
scattering. 

The factorization properties of the scattering amplitude in this domain are conventionally 
described by using the now standard notations of GPDs. For that we rewrite the momenta 
of the particles involved in the process as 

Qi = fii — 2£ n 2 , ki — — rii 

q 2 = --r^-n 1 + (l + On 2 , k 2 = (1 - n 2 , (26) 
(1 + S 

where £ is the skewedness parameter which equals £ = Q\/{2s — Q\) (see also the formula 
(jHJ) which relates s with the total scattering energy W) and the new rii Sudakov light-cone 
vectors are related to the p^s as 

Pi = Y^ n i' P2 = (l + O^a, (27) 

with pi ■ p 2 = rii ■ n 2 = s/2. We also introduce the average "target" momentum P and the 
momentum transfer A 

P = -{q 2 + k 2 ), A = k 2 -q 2 . (28) 

We still restrict our study to the strictly forward case with t = t m i n = — 2£Q%/(1 + £). 

In the region defined by Eq. (|2*ij) we can put Q 2 = inside {...} in the expression in 
Eq. (|13|) . which results, in this approximation, in the formula: 



dzi dz2 <f>(zi)<f>(z 2 ) {....} (29) 


(i-O (1-0 1 1 



dz\ dz 2 4>(zi) 4>(z 2 ) 



^[0 2 (l + O+2^2] ^l[^(l +0+2^2] Z 2 Z! Z X Z 2 



The terms in the {■ • • } are ordered in accordance with diagrams shown in FigO 

Now our aim is to rewrite Eq. (|29| in a form corresponding to the QCD factorization with 

a TDA, as illustrated in FigJTTJ 

For that we look more closely into diagrams contributing to each of the four terms in 

(f29|) from the point of view of such a factorization. For example, the second term in p9"j) . 

corresponding to the second diagram shown in FigJEl suggests the introduction of the new 
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Figure 11: Factorization of the amplitude in terms of a TDA. 



variable x, defined through z 2 = (a? — £)/(! — £), with x G [£, 1], which results in the equality 



dzo 



dx 



21 (x + 



(30) 



i.e. which corresponds to a part of the contribution from the DGLAP integration region of 
the amplitude with a factorized TDA. A similar analysis of all remaining terms permits to 
represent Eq. (|29|) in the form 



1 

/ 



dz\ dz 2 (f>(zi)(f>(z2) {....} 



! 1 

dx J dz\(f){zi) 
l 

©(I > x > 00 



(31) 



1 1 

+ 



Zi(x-0 zi(x + 
x — £ 



e(-e>x>-i) 



1 + x 



with the step fuction 9 (a > x > b) = 9 (a — x)Q(x — b). This factorized expression suggests 
the identification of 



dz\ (f>(zi) 



1 1 

+ 



zt(x-£) z 1 (x + £) j 


as the coefficient function Th (up to a multiplicative factor) , and of 



(32) 



9(l>a;>0 



- 9(-£ >x> -1) 



l + x 



(33) 
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as the Yl ~* Pl TDA. 

To justify this interpretation we start from the hard part T H (zx,x) of the scattering 
amplitude, which appears in FigHU as illustrated at order g 2 in FigJTJJ It equals, for a 





+ 





(X) J»2 fa 

Figure 12: The hard part Th at order g 2 . 




meson built from a quark with a single flavour, 



T H {zi,x) 



■ f 2 r\ ^ I 



ni. 



X 



+ 



_z x {x + £ - it) z 1 (x-£ + ie)_ 



(34) 



and obviously coincides with the hard part of the p— meson electroproduction amplitude. 
The tensorial structure 2£ n 2 + j~+e n i — Pi + Q\l s Vi coincides again with the one present in 
Eq.JBJ). 

Passing to the TDA, let us consider the definition of 72(^2) — ► Pl(^) TDA, T(x, £, t mi „), 
in which we assume that the meson is built from a quark with a single flavour, p 9 L (k 2 ) = qq. 
The vector P = l/2(q 2 + k 2 ) ~ n 2 in our kinematics, and the ray- vector of coordinates is 
oriented along the light-cone vector n = n\j \n\.n 2 ). The non-local correlator defining the 
TDA is given by the formula 

-z/2 

dz~ i t~> \ ~ ieQ " f d y» A *(v) 

' e^(pl(h)\q(-z/2)he ^ q(z/2)\j*(q 2 )) 



2tt 

^ Qq fp 2 

P + Q\ 



e v {q 2 ) (1 + CK + 



Ql 



(35) 



in which we explicitly show the electromagnetic Wilson-line assuring the abelian gauge in- 
variance of the non-local operator. On the contrary, for simplicity of notation, we omit the 
Wilson line required by the non-abelian QCD invariance since it does not play any role in 



this case. Note also that the factor (1 +^)n 2 + 
the tensorial structure of the second term in Eq. 



s(l+0 1 



pu _|_ ¥2.pv corresponds to a part of 
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Now, the simple perturbative calculation of the matrix element in l|3H|l in the lowest order 
in the electromagnetic coupling constant e leads to the expression for T(x, £, t min ) given by 
Eq. (|33|) . The contributing diagrams are drawn in FigIT3l In particular, the contribution to 
the rhs of (|35| proportional to the vector n\ (or p\) corresponds to the contribution coming 
from the expansion of the electromagnetic Wilson line, illustrated by the last diagram in 
FiglU 

Putting all factors together and restoring the flavour structure of the p°, we obtain the 
factorized form involving a TDA, of the expression T a/3 p 2a Pi/3 in Eq. (fT3|) as 

1 1 

T a/ Wi/3 = -if 2 p e\Q 2 u + Q 2 d )g 2 ^ J ' dx J dz, 

-1 o 

(36) 

Note that in this perturbative analysis, only the DGLAP part of the TDA, with 
1 > \x\ > £, contributes. This is a consequence of the support properties of the p— meson 
distribution amplitude. 



+ 



6 Conclusions 

Thus the perturbative analysis of the process 7* 7* — ► p° L p° L in the Born approximation leads 
to two different types of QCD factorization. We have shown that the polarization states of the 
photons dictate either the factorization involving a GDA or involving a TDA. Usually these 
two types of factorizations are applied to two different kinematical regimes. The arbitrariness 
in choosing values of photon virtualities Q\ shows that there may exist an intersection region 
where both types of factorization are simultaneously valid. 

We have restricted our analysis to the case of longitudinally polarized p°— mesons. In this 
way we have avoided the potential problems [10] due to the breaking of QCD factorization 
with GDA or TDA at the end-point region of the distribution amplitudes of transversally 
polarized vector mesons. It would be interesting to find out, whether a factorizing formula 
with GDA or TDA can also be obtained in this case. 

Although we have restricted ourselves, for simplicity, to the forward kinematics, we believe 
that similar results may be obtained in a more general case. Also, we would like to mention 
that our results involving p°— mesons can be generalized to the case of the production of 
a (p + p~)— meson pair. On the theoretical side, in order to preserve the electromagnetic 
gauge invariance of the TDA 7* — > p ± , one should modify the definition of the non-local 
correlator (j35| . This may be done by applying the Mandelstam approach [11], i.e. by 



15 



replacing the electromagnetically gauge invariant corelator (I3H1) by the product of two effective 

oo 

electromagnetically gauge invariant quark fields q(z) exp(i e q J dyn^A^^y)). 

z 
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